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Abstract 

A  review  of  the  experimental  and  computational 
studies  performed  at  NASA  Langley  Research 
Center  (LaRC)  to  support  the  optimization  and 
benchmarking  of  the  hypersonic  aerodynamic 
and  aerothermodynamic  databases  for  the  X-33 
vehicle  is  presented.  A  synoptic  of  the  testing, 
computational,  and  analysis  capabilities  at 
LaRC  applied  to  these  studies  is  given.  Analy¬ 
ses  of  the  hypersonic  aerodynamic  characteris¬ 
tics,  control  surface  effectiveness,  and  reaction 
control  system  effects  are  discussed.  Experi¬ 
mental  measurement  of  the  aerodynamic  heat¬ 
ing  via  the  global  thermographic  phosphor 
technique  and  development  of  a  hypersonic 
boundary-layer  transition  correlation  for  X-33 
is  described.  Computational  results  used  to 
complement  the  experimental  program  and  to 
assess  the  vehicle  aerodynamic  and  aerother¬ 
modynamic  characteristics  in  flight  are  pre¬ 
sented.  The  technical  findings,  impacts,  and 
lessons  learned  from  the  studies  are  discussed. 

1  Overview  of  the  X-33  Program 

1.1  Background 

In  1993,  an  internal  NASA  study  [1]  proposed 
development  of  a  fully  reusable,  rocket  pow¬ 
ered,  single-stage-to-orbit  (SSTO)  vehicle  capa¬ 
ble  of  delivering  25,000  lbs  (including  crew 
members)  to  the  International  Space  Station. 
This  new  launch  system  would  combine  effi¬ 
cient  operations  with  increased  safety  and  pro¬ 
vide  access  to  space  at  greatly  reduced  cost. 
The  “Access-to-Space”  study  identified  critical 
technologies  that  required  development  before  a 
SSTO  reusable  launch  vehicle  (RLV)  could  be 
successfully  flown.  NASA  also  recommended 
that  this  new  system  be  commercially  owned 
and  operated,  unlike  previous  programs.  There¬ 


fore,  an  approach  was  followed  whereby  NASA 
would  work  cooperatively  with  an  industry 
partner  to  mature  the  technology  for  RLV 
through  development  of  a  sub-scale  technology 
demonstration,  i.e.  the  X-33  program. 

The  X-33  program  was  organized  as  an  in¬ 
dustry-led  partnership  with  NASA  in  which  the 
companies  negotiated  development  tasks  with 
the  various  NASA  centers  prior  to  submitting 
proposals.  Industry  would  utilize  government 
facilities  and  share  in  the  program  costs.  A 
Phase  I  competition  for  X-33  was  held  in  1995 
among  the  industry  leaders  (Lockheed-Martin, 
Rockwell,  McDonnell  Douglas/Boeing).  Lock¬ 
heed-Martin  Skunk  Works  (LMSW)  was  ulti¬ 
mately  awarded  the  program  to  continue  into 
Phase  II  and  toward  construction  and  flight  of 
the  X-33  vehicle.  The  X-33  program  (Phase  I 
and  Phase  II)  had  ambitious,  fast-paced  sched¬ 
ules  with  a  total  time  from  development  to  flight 
of  3-4  years.  Development  of  VentureStar  , 
the  full-scale  operational  RLV,  progressed  in 
parallel  with  the  X-33  program  [2]. 


Figure  1.  The  Lockheed-Martin  X-33  and 
early  RLV  compared  to  Shuttle. 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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1.2  Configuration 

The  Lockheed-Martin  concept  for  X-33  (and 
RLV)  follows  a  long  history  of  lifting-body  de¬ 
velopment  [3]  in  the  United  States.  On  X-33, 
canted  fins,  twin  vertical  tails,  and  dual  body 
flaps  are  incorporated  for  aerodynamic  control 
of  the  lifting  body.  The  vehicle  is  powered  by 
two  linear  aerospike  engines  mounted  on  the 
base.  A  rendering  of  X-33  relative  to  the  Space 
Shuttle  and  an  early  RLV  concept  is  presented 
in  Figure  1  while  the  drawings  illustrated  in  Figure 
2  show  more  details  and  give  the  primary  di¬ 
mensions  of  the  vehicle. 
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Figure  2.  Primary  geometric  features  of  the  Lock¬ 
heed-Martin  X-33  vehicle 
(from  Hollis,  et  al.  [4]). 

A  unique  feature  of  the  X-33  is  the  use  of 
metallic  panels  for  the  windward  surface  ther¬ 
mal  protection  system  (TPS).  These  TPS  panels 
vary  in  size  depending  on  body  location  and  are 
arranged  in  a  herringbone  pattern  with  overlap¬ 
ping  edges  to  eliminate  gaps  and  fillers.  A 
sketch  of  the  panels  on  the  windward  surface 
along  with  the  TPS  materials  employed  over  the 
remaining  vehicle  is  given  in  Figure  3.  The  me¬ 
tallic  panels  are  expected  to  provide  benefits  in 
vehicle  operability  relative  to  a  ceramic  TPS 
such  as  employed  on  the  Shuttle  Orbiter. 
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Figure  3.  TPS  material  layout. 


1.3  Flight 

The  X-33  was  designed  as  a  sub-orbital  vehicle 
with  two  basic  flight  trajectories  from  its  launch 
site  at  Edwards  Air  Force  Base  (AFB),  Califor¬ 
nia.  The  first  is  a  high  performance  trajectory 
that  would  reach  250,000  ft  at  Mach  15  and  land 
at  Malmstrom  AFB,  Montanna.  The  second  is  a 
lower  Mach  10  flight  reaching  180,000  ft  to 
Michaels  AFB,  Utah.  Design  of  the  X-33  ther¬ 
mal  protection  system  was  based  on  the  Malm¬ 
strom  trajectory  but  the  first  flights  will  be 
along  the  more  benign  path  to  Michaels.  During 
descent,  the  vehicle  attitude  is  maintained  be¬ 
tween  30  and  40-deg  until  approximately 
Mach  5.  Figure  4  depicts  the  descent  portion  of 
the  two  trajectories  and  shows  the  operating  en¬ 
velopes  of  the  ground-based  facilities  used  in 
the  work  reviewed  herein. 
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Figure  4.  Mach  number  -  Reynolds  number  simula¬ 
tion  capability  (adapted  from  Miller  [5]). 
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The  sub-orbital  flight  demonstrations  of 
X-33  will  test  most  flight  and  operational  as¬ 
pects  of  the  future  RLV;  however,  it  should  be 
noted  the  flights  will  not  provide  aerothermody- 
namic  information  about  the  low  density,  hy¬ 
pervelocity  flight  encountered  during  reentry 
from  low  earth  orbit  (LEO). 

1.4  Role  in  X-33  Program 

The  Aerothermodynamics  Branch  (AB)  at 
NASA  Langley  Research  Center  (LaRC)  has 
played  major  roles  in  Phase  I  and  Phase  II  of  the 
X-33  program.  During  the  1995  Phase  I  study, 
the  AB  simultaneously  supported  the  screening, 
design,  and  optimization  of  three  industry  con¬ 
cepts  over  a  12-month  period.  During  Phase  II, 
formal  Task  Agreements  with  LMSW  were  es¬ 
tablished  whereby  the  AB  would  support  the 
development,  final  design,  and  benchmarking  of 
Lockheed’s  lifting  body  concept.  This  com¬ 
mitment  posed  a  significant  challenge  to  the  AB 
since  resources  were  shared  with  several  com¬ 
peting  programs  (e.g.  X-34,  X-38,  and  X-43). 
Miller  [6,5]  gives  a  general  review  of  the  pro¬ 
grammatic  difficulties  and  the  lessons  learned 
from  the  AB’s  effort  during  this  period  of  in¬ 
tense  work. 

The  purpose  of  the  present  paper  is  to  pro¬ 
vide  a  technical  review  of  several  major  aerody¬ 
namic  and  aerothermodynamic  studies  per¬ 
formed  by  the  AB  in  support  of  the  X-33  devel¬ 
opment.  As  such,  this  report  presents  the  per¬ 
spective  of  NASA  LaRC,  and  the  Aerothermo¬ 
dynamics  Branch  in  particular.  The  technical 
findings,  impact,  and  lessons  learned  from  de¬ 
velopment  of  the  aerodynamic  and  aerothermo¬ 
dynamic  databases  for  the  X-33  vehicle  are  pre¬ 
sented.  A  typical  design  process  supported  by 
the  AB  begins  with  the  screening  of  initial  con¬ 
cepts  and  follows  through  optimization  of  the 
vehicle  shape  and  then  into  benchmarking  of  the 
final  design  as  depicted  in  Figure  5.  Data  from 
these  studies  form  the  basis  for  vehicle  design 
and  flight  and  provide  critical  inputs  to  other 
major  systems  (e.g.  guidance  and  control,  TPS 
design,  structural  loads).  A  final,  but  often- 
unrealized  step,  is  the  validation  of  design  tools 
with  flight-derived  data  as  shown  in  the  figure. 


D»ta 

^ _ Fly  mission  ^ 

Mission  requirements  ] 

_ 4 _ 

i 

1 

Aerothermodynamic  benchmarking 

(  Systems  analysis  ] 

(  \ 

Develop  flight  data  book 

Validate 

Establish  flight  margins 

i 

i 

Wind  j 

-  Aerodynamic/aeroheating 

C  Initial  ^ 

Tunnel  / 

4 

^  concept  ) 

and  / 

C  Outer  mold  ^ 

l CFD/ 

^  lines  (OML)  frozen  J 

i 

- - - 

Aerothermodynamic  screening 

Aerothermodynamic  optimization 

Parametric  ground-based  testing 

Detailed  ground-based  testing 

Engineering  codes 

CFD  codes 

Flyable  throughout  trajectory? 

Optimum  aerodynamic  performance? 

-  Subsonic  - 

hypersonic  aerodynamics 

Yes 

-  Descent,  approach,  landing 

Survivable  ? 

Quantified/optimum  aeroheating 

-Aeroheating 

and  aeroloads? 

Jno 

V 

fno 

Figure  5.  Aerothermodynamic  process  (adapted  from 
Miller  [5]). 

The  next  section  of  this  paper  describes  the 
experimental  and  computational  capabilities  at 
LaRC  that  were  used  to  support  the  design  proc¬ 
ess  for  X-33.  The  remaining  sections  present 
selected  results  that  highlight  the  technical  con¬ 
tributions.  The  paper  focuses  first  on  hyper¬ 
sonic  aerodynamic  testing  (primarily,  pitch 
control  issues)  and  assessment  of  the  reaction 
control  system  (RCS),  followed  by  assessment 
of  global  and  local  aerodynamic  heating,  in¬ 
cluding  a  discussion  of  the  hypersonic  bound- 
ary-layer  transition  correlations  that  were  devel¬ 
oped. 

2  Experimental  Facilities  and  Testing 

2.1  Aerothermodynamic  Facilities  Complex 

The  experiments  reviewed  in  this  paper  were 
conducted  in  four  facilities  of  the  Aerothermo¬ 
dynamic  Facilities  Complex  at  NASA  Langley 
Research  Center:  the  20-Inch  Mach  6  Air  Tun¬ 
nel,  the  31-Inch  Mach  10  Air  Tunnel,  the  20- 
Inch  Mach  6  CF4  Tunnel,  and  the  22-Inch  Mach 
20  Helium  Tunnel.  Table  1  gives  a  summary  of 
the  nominal  reservoir  and  corresponding  free- 
stream  flow  conditions  in  each  facility  along 
with  the  Reynolds  number  range  available.  A 
detailed  description  of  these  facilities  along  with 
their  performance  characteristics,  history,  and 
capabilities  is  given  by  Micol  [7], 
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Table  1.  Nominal  flow  conditions. 


Facility 

Pt,  psi 

T,,°R 

q«*>,  psi 

ReooXlO'6,  ft 1 

Y~ 

Nominal 

Range 

20-Inch  Mach  6  Air 

125 

910 

6 

2.0 

2.0 

0.5  to  8 

1.4 

31 -Inch  Mach  10  Air 

1450 

1800 

10 

2.2 

2.2 

0.5  to  2.1 

1.4 

20-Inch  Mach  6  CF4 

950 

1150 

6 

0.9 

0.4 

0.05  to  0.7 

1.15 

22-Inch  Mach  20  He 

1000 

530 

20 

1.7 

7.5 

2.4  to  22 

1.67 

2.2  Aerodynamic  Testing 

All  of  the  hypersonic  X-33  aerodynamic  tests 
were  performed  with  a  0.007-scale  stainless 
steel  model  designed  and  fabricated  at  LaRC. 
The  model  included  removable  canted  fins, 
body  flaps,  vertical  tails,  and  engine  nozzle 
components.  The  removable  components  were 
interchangeable  with  deflected  control-surface 
parts  or  with  off-blocks  that  were  machined  to 
the  fuselage  contour.  The  model  off-blocks  en¬ 
abled  configuration  buildup  tests  to  investigate 
the  aerodynamics  of  the  baseline  vehicle  (all 
components  installed;  all  deflections  zero)  and 
all  combinations  of  the  fuselage,  control  sur¬ 
faces,  and  engine  nozzle. 

The  six  force  and  moment  components  on 
the  vehicle  were  measured  on  a  top-mounted 
blade  strut  for  angles  of  attack  (a)  over  25  deg 
to  minimize  support  interference  effects.  A 
straight  sting  was  employed  at  lower  angles  of 
attack  with  a  10  deg  overlap  between  the  two 
mounting  systems  to  insure  continuity  of  the 
data.  An  experimental  uncertainty  was  as¬ 
signed  to  all  data  based  on  a  measured  ±0.5% 
uncertainty  in  the  balance  at  full-scale  load. 

2.3  Aerothermodynamic  Testing 

Measurements  to  determine  aerodynamic  heat¬ 
ing  were  performed  using  a  two-color  relative- 
intensity  phosphor  thermography  technique  [8- 
10]  to  enable  nonintrusive  optical  acquisition  of 
the  test  data.  With  this  technique,  silica  ceramic 
wind-tunnel  models  are  slip  casted  [11]  and 
coated  with  a  mixture  of  phosphors  that  fluo¬ 
resce  in  two  regions  of  the  visible  spectrum  (red 
and  green)  when  illuminated  with  ultraviolet 
light.  The  fluorescence  intensity  is  dependent 
on  the  amount  of  incident  light  and  local  tem¬ 
perature  of  the  phosphors.  By  acquiring  fluo¬ 


rescence  intensity  images  of  an  illuminated 
model  exposed  to  the  wind-tunnel  flow,  surface 
temperature  can  be  inferred  for  portions  of  the 
model  in  the  camera’s  view.  A  temperature 
calibration  of  the  system  conducted  prior  to  the 
test  provides  the  data  needed  to  convert  the  two- 
color  images  to  quantitative  values  of  surface 
temperature.  Acquiring  video  images  (tem¬ 
peratures)  at  30  frames/sec  during  the  wind- 
tunnel  run  enables  a  corresponding  global  heat- 
transfer  distribution  to  be  computed.  Compari¬ 
sons  of  heat  transfer  measurements  using  con¬ 
ventional  thin-film  resistance  gauges  and  the 
thermographic  phosphor  technique  have  shown 
[10]  excellent  agreement.  Accuracy  of  the 
phosphor  system  and  data  reduction  is  ±8%  and 
the  overall  accuracy  on  heating  rate  is  ±15%  due 
to  all  factors  [10]. 

This  technique,  which  has  been  widely 
used  at  LaRC,  offers  several  advantages  over 
conventional  test  methods.  Foremost,  the  meas¬ 
urements  provide  a  quantitative  resolution  of 
global  temperatures  and  heating  unlike  discrete 
gauge  measurements.  In  addition,  the  model 
construction,  testing,  and  data  reduction  can  be 
performed  significantly  faster  and  cheaper  than 
other  techniques  allow.  The  global  resolution 
and  rapid  testing  this  method  allows  were  par¬ 
ticularly  important  for  the  X-33  work  in  order  to 
study  transition  fronts  on  a  complex  three- 
dimensional  vehicle  in  a  timely  and  efficient 
manner.  Over  70  ceramic  models  were  con¬ 
structed  to  support  aerothermodynamic  testing 
of  the  X-33  vehicle.  Advances  were  made  in 
fabrication  techniques  of  the  ceramic  models  in 
response  to  program  requirements  (e.g.  models 
with  highly  a  detailed  surface  to  simulate  bow¬ 
ing  of  individual  TPS  panels). 
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3  Computational  Methods 

3.1  Engineering  Code 

Computation  of  surface  heating,  surface  stream¬ 
lines,  and  boundary  layer  parameters  was  per¬ 
formed  using  the  LATCH  (Langley  Approxi¬ 
mate  Three-Dimensional  Convective  Heating) 
code  [12].  LATCH  is  based  on  the  axisymmet- 
ric  analog  and  uses  approximate  integral  meth¬ 
ods  for  general  three-dimensional  boundary 
layer  results.  The  code  is  coupled  with  inviscid 
(Euler)  solutions  to  provide  the  boundary  layer 
edge  conditions.  LATCH  employs  an  engi¬ 
neering  technique  [13]  to  calculate  heating  on 
three-dimensional  hypersonic  vehicles  and  has 
been  shown  [12]  to  be  in  good  agreement  with 
both  experimental  and  Navier-Stokes  results. 

3.2  Euler  Code 

The  FELISA  code  [14,15]  combines  a  series  of 
programs  that  generate  unstructured  tetrahedral 
grids  over  three-dimensional  vehicles  and  solve 
the  steady  Euler  equations.  This  proves  to  be  an 
invaluable  combination  since  it  enables  very 
rapid  modeling  and  solutions  for  complex  con¬ 
figurations.  Algorithms  exist  in  FELISA  to 
solve  flows  ranging  from  transonic  to  hyper¬ 
sonic  speeds  with  gas  chemistry  models  for  air 
(perfect  gas  and  equilibrium),  CF4,  and  the  Mars 
atmosphere  (CO2). 

3.3  Navier-Stokes  Codes 

Two  Navier-Stokes  codes  were  routinely  ap¬ 
plied  in  work  reviewed  herein.  The  LAURA 
(Langley  Aerothermodynamic  Upwind  Relaxa¬ 
tion  Algorithm)  code  [16]  uses  a  finite-volume 
shock-capturing  approach  to  solve  high-speed 
(Moo>3)  steady  viscous  and  inviscid  flow  prob¬ 
lems.  The  algorithm  incorporates  a  point- 
implicit  relaxation  scheme  to  obtain  solutions 
efficiently  on  multi-processor  and  massively 
parallel  computers.  The  code  has  been  success¬ 
fully  applied  to  a  full  range  of  hypersonic  vehi¬ 
cles  and  flight  conditions  during  its  evolution 
over  the  past  decade.  LAURA  includes  models 
for  perfect  gas,  equilibrium  air,  and  thermal  and 
chemical  nonequilibrium  air  in  addition  to  mod¬ 


els  for  CF4  and  a  Mars  atmosphere.  A  unique 
feature  of  the  LAURA  code  is  the  ability  to  per¬ 
form  one-dimensional  grid  adaption  in  order  to 
resolve  high  gradients  in  the  boundary  layer  and 
across  a  bow  shock. 

Like  LAURA,  the  GASP  (General  Aero¬ 
dynamic  Simulation  Program)  code  [17]  em¬ 
ploys  a  finite- volume,  shock-capturing  approach 
to  solve  both  inviscid  and  viscous  flow  fields  on 
structured  grids.  GASP  employs  a  wide  range 
of  algorithms  and  models  to  solve  flows  from 
subsonic  through  hypersonic  conditions. 
Equivalent  thermochemical  models  exist  in  the 
GASP  code  in  comparison  with  LAURA.  The 
methodology  for  one-dimensional  grid  adaption 
employed  by  LAURA  has  been  incorporated 
into  the  GASP  solution  procedure. 

3.4  Grid  Topology  and  Generation 

The  complex  geometry  in  the  base  region  of 
X-33  posed  a  challenge  for  grid  generation  and 
analysis  with  detailed  CFD  (Navier-Stokes  and 
Euler).  Consequently,  two  models  of  the  base 
region  emerged;  the  first  attempted  to  capture 
the  full  geometric  complexity  of  the  vehicle  and 
the  second  used  an  extended  fuselage  to  sim¬ 
plify  the  problem.  Figure  6  illustrates  the  fuse¬ 
lage  extension  that  covered  the  base  region  and 
permitted  solutions  over  the  entire  surface  (in¬ 
cluding  body  flap  and  canted  fin)  without  cal¬ 
culation  of  the  complicated  wake  flow.  A  lim¬ 
ited  number  of  computations  were  performed 
without  this  fuselage  extension  but  the  majority 
of  cases  used  the  simplified  geometry. 


Figure  6.  Computational  surface  model  of  X-33  ge¬ 
ometry. 
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Structured  grid  generation  has  traditionally 
been  a  difficult  and  time-consuming  process  in 
the  application  of  Navier-Stokes  codes.  How¬ 
ever,  the  arduous  requirements  for  grid  genera¬ 
tion  on  X-33  promoted  improvements  in  the 
grid  generation  process;  in  particular,  the  capa¬ 
bility  to  smooth  and  manipulate  complex  struc¬ 
tured  grids.  Alter  [18,19]  refined  the  methodol¬ 
ogy  and  software  that  performed  grid  smoothing 
and  manipulation  into  the  VGM  (Volume  Grid 
Manipulator)  code  [20]  during  this  period.  The 
VGM  code  evolved  from  a  simple  grid  conver¬ 
sion  tool  to  a  sophisticated  scripting  program 
that  proved  invaluable  for  grid  quality  im¬ 
provements,  grid  repairs,  and  batch  processing 
of  grids  among  other  features. 

4  Contributions  to  the  X-33  Program 

4.1  Aerodynamics 

Aerodynamic  screening  of  the  Lockheed  X-33 
configuration  was  completed  and  iterations  to 
optimize  the  aerodynamic  stability  and  control 
were  initiated  during  the  Phase  I  competition. 
Phase  II  work  was  mainly  intended  to  generate 
the  benchmark  data  to  be  included  in  the  vehicle 
aerodynamic  flight  databook.  The  focus  of  this 
review  concerns  the  testing  performed  in  the 
LaRC  hypersonic  facilities  during  Phase  II. 
These  tests  represented  a  major  contribution  to 
the  overall  experimental  aerodynamic  database 
achieved  via  tests  at  other  Langley  facilities 
(e.g.  the  Low  Pressure  Turbulence  Tunnel 
(LTPT),  the  14x22-Foot  Wind  Tunnel,  the  Uni¬ 
tary  Plan  Wind  Tunnel  (UPWT),  and  16-Foot 
Transonic  Tunnel)  and  wind  tunnels  across  the 
United  States  to  cover  the  Mach  number  and 
Reynolds  number  ranges  of  flight.  The  applica¬ 
tion  of  CFD  in  these  efforts  was  used  to  support 
the  wind  tunnel  testing  through  comparison  with 
the  measured  data  and  by  interrogating  the  so¬ 
lutions  to  better  understand  the  fluid  dynamics. 
CFD  applications  to  the  flight  vehicle  were  also 
made  to  assess  extrapolation  of  the  wind  tunnel 
data  to  flight  conditions. 

During  the  series  of  hypersonic  tests,  over 
600  runs  were  made  in  the  four  LaRC  aerother- 
modynamic  facilities.  Tests  were  performed  for 


angles  of  attack  from  -4  to  50  deg  and  sideslip 
angles  from  -4  to  4  deg  to  measure  the  aerody¬ 
namic  forces  and  moments.  The  effects  of  con¬ 
trol  surface  deflections  were  studied  by  testing 
body  flap  settings  between  -15  and  30  deg  and 
elevon  settings  between  -30  to  30  deg.  No  rud¬ 
der  deflections  were  considered  at  the  hyper¬ 
sonic  conditions  studied.  By  conducting  tests  in 
both  Mach  6  and  Mach  10  air  flows  a  Mach 
number  effect  could  be  inferred  in  the  aerody¬ 
namic  database.  Similarly,  tests  at  Mach  6  in 
both  air  and  CF4  were  used  to  simulate  real  gas 
effects  through  the  variation  in  shock  density 
ratio  produced  in  the  different  test  gases.  Tests 
showed  a  measurable  effect  over  the  Reynolds 
number  range  in  each  facility,  but  the  effect  was 
quite  small  and  is  not  discussed  in  this  paper. 
Tests  in  Mach  20  helium  were  used  for  screen¬ 
ing  and  rapid  assessment  with  stereo¬ 
lithography  (SLA)  models  and  were  not  part  of 
the  aerodynamic  database.  CFD  solutions  using 
the  LAURA  and  GASP  Navier-Stokes  codes 
and  the  FELISA  Euler  code  were  obtained 
[21,22]  over  the  high  angle  of  attack  range  for 
both  Mach  6  air  and  CF4  and  Mach  10  air.  Vis¬ 
cous  solutions  with  deflected  body  flap  controls 
were  also  modeled  in  the  computations. 

A  complete  review  of  the  hypersonic  aero¬ 
dynamic  characteristics  established  in  the  LaRC 
tests  has  been  presented  by  Murphy,  et  al.  [23]. 
In  general,  the  tests  showed  the  vehicle  to  have 
a  linear  lift  curve  for  a<40  deg  and  a  hypersonic 
L/D  of  approximately  1  at  trim  condition.  A 
good  comparison  between  measurement  and 
computation  of  axial  force  (Ca)  and  normal 
force  (Cn)  was  found  across  the  angle  of  attack 
and  Mach  number  range.  The  experiments 
showed  the  vehicle  to  be  roll  stable  for  a  >  4 
deg  and  to  posses  positive  dihedral  effect.  Not 
surprisingly,  the  vehicle  displayed  directional 
instability  for  all  angles  of  attack  at  these  hyper¬ 
sonic  conditions  due  to  the  high  incidence  and 
shadowing  of  the  vertical  tail  surfaces.  This  is  a 
common  trait  of  a  winged  or  lifting-body  entry 
vehicle  that  is  typically  mitigated  with  reaction 
control  jets. 
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4.1.1  Longitudinal  Control 
Hypersonic  pitching  moment  (Cm)  is  presented 
in  Figure  7  as  a  function  of  angle  of  attack  for 
Mach  6  and  10  conditions  in  air.  As  shown  in 
the  figure,  the  vehicle  exhibits  longitudinal  in¬ 
stability  for  a  <  10°,  neutral  stability  for  10°  <  a 
<  24°,  and  then  marginal  stability  for  higher  an¬ 
gles  of  attack.  The  measured  data  and  the  CFD 
predictions  display  a  crossover  of  pitching  mo¬ 
ment  between  the  Mach  6  and  Mach  10  condi¬ 
tions  where  the  crossover  was  measured  to  oc¬ 
cur  around  a=30°  but  is  predicted  by  CFD  to 
occur  around  a=42°.  A  parametric  CFD  study 
[21]  confirmed  a  consistent  trend  of  decreasing 
stability  (and  crossover)  with  decreasing  Mach 
number  for  the  X-33,  and  pitching  moment  data 
from  other  configurations  (i.e.  winged  bodies 
and  lifting  bodies)  tested  in  the  LaRC  Mach  6 
and  10  air  tunnels  show  similar  Mach  number 
effects.  Subsequent  configuration  build-up 
studies  on  X-33  identified  the  fuselage  alone  as 
the  source  of  the  Mach  number  effect,  although 
the  exact  cause  has  not  been  determined. 


The  largest  disagreement  between  meas¬ 
ured  and  predicted  pitching  moment  occurs  in 
the  angle  of  attack  range  above  a=30°,  particu¬ 
larly  for  Mach  6  where  the  prediction  is  outside 
the  experimental  uncertainty  band.  The  dis¬ 
crepancy  for  Mach  6  naturally  raised  concerns 
with  both  the  CFD  predictions  and  wind  tunnel 
measurements  since  the  baseline  vehicle  (zero 
control  deflections)  was  less  stable  with  de¬ 
creasing  Mach  number.  It  was  thought  that  the 
failure  to  model  the  base  and  wake  region  in  the 
computations  could  present  some  error;  how¬ 
ever,  a  CFD  solution  with  the  full  laminar  wake 
at  a=36  deg  showed  little  effect  of  the  base  flow 
on  the  baseline  pitching  moment  (Figure  7). 
Comparing  results  from  the  configuration  build¬ 
up  tests  with  the  moments  computed  over  indi¬ 
vidual  components  isolated  the  source  of  the 
disagreement  to  the  fuselage  alone. 

One  additional  study  was  performed  using 
a  sparsely  instrumented  pressure  model  in  a 
rapid  attempt  to  understand  the  differences  in 
the  experimental  and  computational  results  for 
Mach  6  and  10  air.  The  pressure  tests  suggested 
that  differences  existed  in  the  measured  and 
predicted  pressure  on  the  wind  side  forebody  in 
an  area  of  overexpansion  and  recompression 
downstream  of  the  nose  cap.  However,  a  con¬ 
sistent  explanation  for  both  Mach  6  and  10 
could  not  be  drawn  and  the  results  were  incon¬ 
clusive.  A  high  fidelity,  highly  instrumented 
pressure  model  is  required  to  provide  the  neces¬ 
sary  information  to  explain  the  fluid  dynamic 
phenomena  causing  the  Mach  number  effect  and 
resolve  the  disagreement  between  experiment 
and  computation.  Nevertheless,  it  is  important 
to  note  that  the  discrepancies  in  pitching  mo¬ 
ment  coefficient  due  to  variation  in  Mach  num¬ 
ber  from  6  to  10  are  small  (on  the  order  of 
0.008)  and  that  the  X-33  body  flaps  produce 
sufficient  authority  to  control  the  vehicle  as  dis¬ 
cussed  in  the  next  section. 


Figure  7.  Comparison  of  measured  and  predicted 
pitching  moment  (baseline),  (adapted  from  Murphy, 
etal.  [23]). 
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4.1.2  Control  Surface  Effectiveness 
Rear-mounted  body  flaps  (Figure  2)  provide  the 
primary  pitch  control  on  the  X-33  vehicle  dur¬ 
ing  unpowered  descent.  The  effectiveness  of 
these  flaps  is  illustrated  in  Figure  8,  which  shows 
the  pitching  moment  curves  for  positive  flap 
deflections  (8bf)  of  10  and  20  deg  at  Mach  10. 
Mach  6  measurements  (not  shown)  are  similar 
but  slightly  less  stable  and  slightly  under  pre¬ 
dicted  by  the  CFD  (as  was  observed  for  the 
baseline  case). 


Figure  8.  Comparison  of  measured  and  predicted 
pitching  moment  with  deflected  body  flaps  (from 
Hollis,  et  al.  [4]). 

From  a  programmatic  viewpoint,  the 
curves  show  that  the  vehicle  has  more  than  suf¬ 
ficient  control  authority  for  stability  and  trim 
and  the  previously  discussed  uncertainties  in 
pitching  moments  translate  to  approximately  2-3 
deg  of  body  flap  deflection.  From  a  fluid  dy¬ 
namic  perspective,  the  curves  portray  degrada¬ 
tion  in  flap  effectiveness  as  angle  of  attack 
and/or  body  flap  deflection  are  increased.  This 
degradation  is  due  primarily  to  an  interaction 
between  the  bow  shock  and  deflected  flap  shock 
that  produces  a  reflected  expansion  fan  over  the 
flap  surface.  The  CFD  predictions  captured  this 
interaction  and  the  corresponding  influence  on 
aerodynamic  characteristics  very  well,  thereby 


complementing  the  experimental  data  and  pro¬ 
viding  invaluable  information  toward  under¬ 
standing  the  structure  of  this  complex  flow  field. 

4.1.3  Reaction  Control  System  (RCS) 

Reaction  control  jets  are  an  integral  component 
of  space  transportation  systems  since  they  are 
required  for  on-orbit  maneuver  and  provide 
control  when  aerodynamic  forces  are  insuffi¬ 
cient  during  periods  of  low-density  flight  or 
high  vehicle  incidence.  In  the  case  of  X-33,  the 
planned  trajectory  remains  within  the  sensible 
atmosphere  but  includes  high  angle  of  attack 
bank  maneuvers  that  require  a  RCS.  A  primary 
issue  in  RCS  design  is  the  possible  aerodynamic 
influence  of  interactions  produced  between  a 
firing  jet  plume,  the  surrounding  flow  field,  and 
the  vehicle  surface.  These  RCS  interactions  are 
difficult  to  accurately  predict  at  flight  conditions 
and  can  be  significant  as  was  learned  for  the 
Shuttle  Orbiter  during  the  first  flight  (STS-1) 

[24] , 

A  system  of  eight  aft-end  thrusters  (4  per 
side)  was  designed  for  the  X-33  vehicle  to  pro¬ 
vide  the  needed  pitch,  yaw,  and  roll  control  as 
shown  in  Figure  9.  A  series  of  experimental  tests 
was  performed  in  the  Mach  10  air  and  Mach  6 
air  tunnels  and  the  LaRC  UPWT  at  Mach  4  to 
determine  the  level  of  interaction  effects  on  the 
X-33  aerodynamics.  Tests  were  performed  with 
a  1%  scale  model  and  unheated,  high-pressure 
nitrogen  gas  was  used  to  simulate  the  thruster’s 
jet  plumes.  Chamber  pressure  to  the  thruster 
nozzles  was  varied  from  400  to  1200  psi  to 
simulate  the  change  in  jet  plume  shape  with  al¬ 
titude  during  flight.  Three  sets  of  conical  noz¬ 
zles  with  exit  half-angles  between  10  and  30  deg 
were  tested  to  simulate  the  change  in  plume 
shape  during  flight  due  to  the  different  specific 
heat  ratios  (y)  in  the  actual  jet  and  external  flow 
field.  Variation  of  these  simulation  parameters 
was  guided  by  Shuttle  flight  data  where  Scallion 

[25]  showed  matching  momentum  ratio  of  the 
jet  and  freestream  gave  the  best  correlation  be¬ 
tween  wind  tunnel  data  and  flight.  RCS  inter¬ 
action  on  deflected  control  surfaces  (body  flap, 
elevon,  and  rudder)  was  also  examined  during 
the  experiment.  Results  from  these  tests  pro¬ 
vided  the  vehicle  designers  with  a  complete  set 
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of  force  and  moment  increments  on  the  baseline 
aerodynamics  that  could  be  used  to  evaluate  in¬ 
teraction  effects  and  design  control  systems.  It 
was  found  that  side-firing  jets  located  below  the 
X-33  canted  fins  (jet  nozzle  number  2  in  Figure 
9)  produced  uncontrollable  rolling  moments  due 
to  interactions  with  the  fin  and  these  yaw- 
control  jets  were  subsequently  moved  to  the  ve¬ 
hicle  base. 


Figure  9.  Sketch  of  wind  tunnel  model  for  reaction 
control  system  (RCS)  tests. 

4.2  Aerothermodynamics 

The  contributions  to  the  aerothermodynamic 
design  for  X-33  summarized  in  this  paper  in¬ 
clude  continued  development  and  validation  of 
the  tools  used  to  predict  the  aeroheating  envi¬ 
ronment  and  the  application  of  these  tools  for 
design  and  benchmark  data.  Development  of 
aerothermodynamic  tools  progressed  on  two 
fronts.  Firstly,  use  of  thermographic  phosphors 
for  heating  measurements  eliminated  the  usual 
situation  where  experimental  aerothermody¬ 
namic  analysis  lags  the  aerodynamic  analysis 
due  to  the  complexity  of  aerothermodynamic 
models  and  instrumentation.  The  use  of  ceramic 
phosphor  models  provided  an  opportunity  to 
conduct  both  type  of  analyses  in  parallel.  This 
is  an  important  step  since  aerodynamic  optimi¬ 
zations  often  change  the  outer  mold  line  shape 
that  in  turn  requires  more  aeroheating  assess¬ 
ment.  The  second  area  of  development  per¬ 
tained  to  the  prediction  of  hypersonic  boundary 
layer  transition  on  X-33.  Although  progress 


has  been  made  in  predicting  transition  from  the¬ 
ory,  in  practice  vehicle  designers  still  rely  on 
semi-empirical  techniques  such  as  used  for 
Shuttle.  However,  advances  in  CFD  tools  and 
global  experimental  techniques  (e.g.  thermo¬ 
graphic  phosphors)  have  significantly  improved 
these  techniques  both  in  accuracy  and  in  time 
required  for  analysis. 

Over  1100  runs  were  performed  in  the 
LaRC  AFC  during  the  coarse  of  this  work  to 
provide  essentially  all  of  the  experimental  aero¬ 
heating  data  for  the  X-33  program.  Tests  were 
performed  [26,27]  in  the  Mach  6  air  facility  at 
angles  of  attack  from  -5  to  40  deg  to  simulate 
both  ascent  and  descent  conditions.  Body  flap 
deflections  between  0  and  20  deg  were  studied 
and  freestream  Reynolds  numbers  based  on 
body  length  were  varied  from  0.8  to  6.6  million 
in  the  facility.  Thermographic  phosphor  meas¬ 
urements  were  used  to  capture  heating  distribu¬ 
tions  for  the  vehicle  acreage  as  well  as  for  spe¬ 
cialized  areas  of  heating.  The  following  sec¬ 
tions  describe  the  results  of  these  global  and  lo¬ 
cal  aeroheating  studies  including  the  supporting 
computations. 

4.2.1  Global  Aeroheating 
Aeroheating  analyses  were  performed  [28]  dur¬ 
ing  X-33  Phase  I  using  a  limited  number  of 
CFD  solutions  at  points  along  the  flight  trajec¬ 
tory  coupled  with  approximate  methods  to 
populate  an  aeroheating  database.  A  similar 
method  was  employed  [29]  in  Phase  II  to  refine 
the  aeroheating  LMSW  X-33  database.  Valida¬ 
tion  of  the  CFD  codes  used  to  anchor  the  data¬ 
bases  (i.e.  LATCH  and  GASP)  was  accom¬ 
plished  by  comparison  of  computed  heating 
rates  to  those  measured  in  the  hypersonic  wind 
tunnel  tests.  Additional  comparison  of  results 
from  the  different  computational  tools  (i.e. 
code-to-code  comparisons)  at  both  wind  tunnel 
and  flight  conditions  added  confidence  to  the 
predictions  in  both  environments.  Ultimately, 
comparison  of  predicted  results  (whether  ex¬ 
perimental  or  computational)  with  measured 
flight  data  will  provide  the  final  verification. 
Prior  to  obtaining  such  flight-test  data  for  the 
X-33,  comparisons  with  existing  flight  data  (e.g. 
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Shuttle  orbiter)  is  used  as  documented  for  the 
LAURA  code  in  [30-32], 

Figure  10  illustrates  a  typical  comparison 
between  experimental  and  computational  heat¬ 
ing  rate  along  the  X-33  windward  centerline.  In 
this  case,  both  laminar  and  turbulent  GASP  so¬ 
lutions  are  compared  with  thermographic  phos¬ 
phor  data  at  Mach  6  for  various  values  of  Rey¬ 
nolds  number.  The  data  exhibit  a  transitional 
behavior  beginning  at  the  aft  end  and  moving 
forward  as  Reynolds  number  is  increased.  The 
computed  heating  rates  are  shown  to  agree  well 
with  the  laminar  experimental  data  and  with  the 
turbulent  data  obtained  by  tripping  the  flow. 


Figure  10.  Comparison  of  measured  and  predicted 
heating  along  windward  centerline  (from  Hollis,  et  al. 

[4]). 

Figure  11  compares  measured  lateral  distri¬ 
butions  of  windward  heating  data  with  predic¬ 
tions  using  the  GASP  and  LATCH  codes  at  two 
axial  locations.  Agreement  is  good  with  the  ex¬ 
ception  that  LATCH  over  predicted  heating  in 
the  centerline  region  as  shown  in  Figure  11(a). 
LATCH  predictions  in  this  region  compared 
better  with  GASP  predictions  and  measurement 
as  angle  of  attack  was  increased.  For  the  most 
rearward  axial  station  (Figure  11(b)),  the  predic¬ 
tions  and  experimental  data  are  in  reasonably 


good  agreement  across  the  span.  Note  that  the 
lateral  distribution  at  this  location  includes  both 
the  windward  fuselage  and  the  canted  fin  sur¬ 
faces.  LATCH  does  a  remarkable  job  capturing 
this  complex  heating  pattern. 


y/L 

b)  x/L  =  0.90 


Figure  11.  Comparison  of  measured  and  predicted 
heating  along  lateral  windward  lines  (from  Hollis,  et 
al.  [21]) 

In  other  comparisons,  it  was  found  that 
laminar  CFD  computations  generally  predicted 
the  leeside  measured  heating  rate  within  ex¬ 
perimental  uncertainty.  Finally,  a  typical  com¬ 
parison  between  computations  via  the  GASP 
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and  LAURA  codes  at  flight  conditions  is  shown 
in  Figure  12  where  contours  of  radiation  equilib¬ 
rium  temperatures  show  excellent  agreement. 
(The  reader  is  referred  to  [33,21,4]  for  addi¬ 
tional  examples  of  data  and  code-to-code  com¬ 
parisons). 

Mach  =  8.83 


Figure  12.  Comparison  of  predicted  heating  for  flight 
(from  Hollis,  et  al.  [21  ]). 

A  successful  effort  to  extrapolate  thermo¬ 
graphic  phosphor  heating  measurements  directly 
from  the  wind  tunnel  (i.e.  immediately  follow¬ 
ing  a  tunnel  run)  to  a  full-scale  vehicle  at  flight 
conditions  was  demonstrated  for  the  X-34 
[10,34],  At  present,  this  technique  is  capable  of 
providing  radiation  equilibrium  wall  tempera¬ 
tures  in  flight  for  perfect  gas,  laminar  or  turbu¬ 
lent  flows  as  part  of  the  wind  tunnel  data  reduc¬ 
tion  process.  An  example  of  this  technique  ap¬ 
plied  to  X-33  is  presented  in  Figure  13  for  a 
40-deg  angle  of  attack  condition.  Here  a  Mach 
6  wind  tunnel  measurement  at  length  Reynolds 
number  of  Re™  l=2x106  was  extrapolated  to  a 
flight  environment  at  Mach  6.6,  Re™,L  =5xl06, 
and  146730  ft  (44.7  km)  altitude.  The  figure 
compares  the  surface  temperatures  extrapolated 
along  the  windward  centerline  with  predicted 
results  from  GASP  at  the  flight  condition. 
Temperatures  from  the  two  predictions  are  gen¬ 
erally  within  50  °F.  Obtaining  such  information 
during  the  rapid,  global  aeroheating  tests  af¬ 
forded  by  the  phosphor  technique  is  extremely 
useful  to  vehicle  designers.  Work  to  further 
validate  and  extend  the  extrapolation  procedure 


to  real  gas  flows  is  ongoing  with  the  goal  to  de¬ 
velop  a  rapid  technique  for  TPS  material  selec¬ 
tion,  split  line  definition,  and  material  sizing. 
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Figure  13.  Extrapolation  of  laminar  experimental 
heating  data  to  flight  condition  (from  Horvath,  et  al. 
[26J). 


4.2.2  Local  Aeroheating 
Three  areas  of  prominent  heating  on  X-33  war¬ 
ranted  detailed  investigation;  the  deflected  body 
flap  surfaces,  the  canted  fin  surface  and  leading 
edge,  and  the  aerospike  nozzle.  Results  of  those 
investigations  [26,28]  are  briefly  summarized 
here.  Horvath,  et  al.  [26]  used  phosphor  ther¬ 
mography,  schlieren  images,  and  oil  flow  pat¬ 
terns  to  investigate  the  flow  separation  and  re¬ 
attachment  around  the  body-flap  hinge  line  and 
the  interaction  effects  on  body  flaps  with  10  and 
20-deg  deflections.  Ratios  of  deflected  to  unde¬ 
flected  flap  heating  were  derived  for  use  in  es¬ 
timating  heating  augmentation  over  the  unde¬ 
flected  level.  The  augmentation  at  angles  of 
attack  between  30  and  40  deg  was  found  to  in¬ 
crease  linearly  with  deflection  angle,  with 
maximum  values  of  5  to  7  over  a  laminar  unde¬ 
flected  reference  and  2  to  3  over  a  turbulent  ref¬ 
erence.  The  augmentation  factors  were  found 
insensitive  to  the  state  of  the  approaching 
boundary  layer  (laminar  or  turbulent),  which 
suggests  that  the  flow  was  transitional,  if  not 
fully  turbulent,  over  the  deflected  surfaces. 
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fan  impingement 


a)  Computed  heating  contours,  Reoo,L  = 
3.3x10*  (from  Hollis,  et  al.  [4]) 


b)  Oil  flow  image,  Reoo,L  =  1.7xl06 
(from  Horvath,  et  al.  [26]) 


Figure  14.  Shock  interaction  effects  on  deflected  body 
flap,  M«,=6,  tx=40  deg. 

The  outboard  location  of  the  X-33  body 
flaps  results  in  the  deflected  flap  shock  being  in 
close  proximity  to  the  windward  bow  shock. 
The  interaction  of  these  shocks  produces  a  re¬ 
flected  expansion  fan  that  was  shown  (Figure  8) 
to  cause  degradation  of  the  flap  aerodynamic 
effectiveness.  It  is  not  surprising  that  the  inter¬ 
action  adversely  affects  the  aerodynamic  heat¬ 
ing  over  the  flap  surface  as  well.  Inspection  of 
the  phosphor  temperature  images  reveals  a  local 
heating  peak  produced  by  a  thinning  boundary 
layer  where  the  expansion  fan  impacts  the  sur¬ 


face.  Solutions  obtained  with  the  GASP  code 
(Figure  14(a))  depict  this  local  peak  heating  in 
addition  to  the  separation  and  reattachment 
around  the  hinge  line.  These  predicted  heating 
contours  were  in  qualitative  agreement  with  the 
phosphor  images  although  the  experiment  sug¬ 
gested  the  flow  downstream  of  reattachment 
was  transitional/turbulent  and  calculations  as¬ 
sumed  a  laminar  flow.  Oil  flow  images  such  as 
shown  in  Figure  14(b)  clearly  reveal  the  effects  of 
body  flap  impingement  and  show  good  correla¬ 
tion  with  the  heating  contours.  (The  interaction 
effect  on  the  body  flaps  was  not  observed  in 
studies  on  earlier  X-33  configurations  with 
shorter  length  flaps). 

A  second  area  of  shock  interaction  exists 
on  X-33  between  the  bow  shock  and  the  shock 
formed  on  the  leading  edge  of  the  canted  fin. 
Horvath,  et  al.  [26]  found  that  this  interaction 
occurred  at  all  angles  of  attack  tested  and  most 
closely  resembled  a  Type  VI.  Although  heating 
to  the  fin  leading  edge  in  the  region  of  the  shock 
interaction  could  not  be  discerned  from  the 
phosphor  measurements,  the  global  images  over 
the  fin  surface  did  offer  useful  information.  At 
low  angles  of  attack,  the  shock  interaction  pro¬ 
duced  local  elevated  heating  rates  on  the  upper 
fin  surfaces  due  to  the  negative  incidence  of  the 
fin.  This  region  of  elevated  heating  moved  to 
the  lower  fin  surface  with  increasing  angle  of 
attack  but  the  corresponding  level  was  also  di¬ 
minished.  The  disturbed  flow  downstream  of 
the  shock  interaction  was  sensitive  to  Reynolds 
number  as  expected,  with  transition  over  the  fin 
surface  observed  for  Re™, L  >  3.3  x  106  at  all  an¬ 
gles  of  attack. 

The  aerospike  engine  nozzles  on  X-33  ex¬ 
tend  well  into  the  wake  flowfield  (see  Figure  2) 
and  are  unshielded  from  impingement  of  sepa¬ 
rated  flow  off  the  windward  fuselage.  An  as¬ 
sessment  of  aeroheating  to  the  nozzles  and  base 
was  critical  given  the  high  incidence  angles 
during  descent  and  both  experimental  [26]  and 
computational  studies  [28]  were  undertaken  for 
this  purpose.  Data  from  thermographic  phos¬ 
phor  measurements  indicated  that  heating  rates 
on  the  nozzle  reached  10%  of  the  nose  stagna¬ 
tion  heating  in  the  area  of  impingement.  In 
comparison,  calculations  at  the  flight  peak 
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heating  point  using  the  LAURA  code  predicted 
heating  rates  on  the  order  of  5%  of  the  nose 
stagnation  value  using  a  radiation  equilibrium 
wall  temperature  and  a  surface  emissivity  of 
0.15  on  the  nozzle.  The  complexity  of  the  CFD 
analysis  needed  to  obtain  these  results  is  illus¬ 
trated  in  Figure  15  which  shows  that  250  Cray  C- 
90  hours  were  required  to  obtain  a  solution  and 
that  65%  of  the  total  time  was  required  to  re¬ 
solve  the  wake.  This  was  a  critical  calculation 
that  led  to  the  design  of  a  system  that  circulates 
hydrogen  fuel  through  the  structure  during  de¬ 
scent  to  actively  cool  the  nozzle. 


250  hours  for  solution 

a)  Computer  resources  required 


=11.4 

a  =  34.8  deg 
Re^  L  =  3  x  1 0® 


b)  Computed  streamlines  and  pressure  con¬ 
tours 

Figure  15.  Computation  of  wake  flow  and  impinge¬ 
ment  on  aerospike  nozzle  (from  Miller  [6J). 

4.2.3  Boundary  Layer  Transition 
A  majority  of  the  aerothermodynamic  wind  tun¬ 
nel  testing  performed  for  X-33  and  a  significant 
computational  effort  was  directed  toward  under¬ 
standing  and  developing  a  correlation  for  hyper¬ 
sonic  boundary  layer  transition  on  the  windward 
fuselage.  The  steps  taken  to  investigate  wind 
side  transition  were  straightforward  but  the  or¬ 
der  in  which  the  results  were  applied  was 
somewhat  reversed  due  to  the  fast-paced  nature 
of  the  program.  To  explain  further,  the  vehicle 
designers  used  a  “worst  case”  aeroheating  tra¬ 


jectory  to  design  the  TPS  prior  to  finalizing  the 
transition  correlation  by  assuming  that  turbu¬ 
lence  would  occur  below  a  certain  point  during 
descent.  Once  a  correlation  was  developed, 
then  the  original  assumption  regarding  transition 
onset  would  be  revised,  if  necessary,  and  any 
roughness  constraints  applied.  It  was  assumed 
in  the  program  that  the  flow  adjacent  to  the  lee- 
side  and  body  flap  surfaces  would  always  be 
turbulent  while  the  flow  over  the  canted  fin 
would  transition  30  seconds  prior  to  the  flow 
over  the  fuselage.  As  a  result,  the  work  re¬ 
viewed  in  this  section  focused  exclusively  on 
fuselage  transition. 

The  first  step  in  developing  a  transition 
criterion  was  to  perform  thermographic  phos¬ 
phor  tests  on  a  smooth  model  in  the  20-Inch 
Mach  6  tunnel  over  a  range  of  Reynolds  num¬ 
bers  and  attitudes  and  attempt  to  correlate  the 
transition  fronts  with  the  parameter  Ree/Me. 
Here,  Reeis  the  momentum  thickness  Reynolds 
number  and  Me  the  boundary  layer  edge  Mach 
number,  both  of  which  were  extracted  from 
LATCH  solutions.  As  an  example.  Figure  16(a) 
presents  a  phosphor  image  of  the  transition  front 
on  the  windward  fuselage  at  40-deg  angle  of 
attack  while  Figure  16(b)  shows  a  correlation  of 
this  front  for  a  value  of  Ree/Me  =  300.  Applying 
this  process  to  all  test  conditions  yielded  values 
between  250  and  325  and  an  average  value  of 
Ree/Me  =  285  was  chosen  to  represent  the  whole 
dataset  thereby  giving  a  constant  “smooth  body” 
ratio. 


a)  Thermographic  phosphor  image 
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Figure  16.  Correlation  of  boundary  layer  transition 
on  “smooth”  forebody  (from  Thompson,  et  al.  [35]). 

The  smooth-body  transition  tests  provided 
an  upper  bound  on  the  Ree/Me  parameter  since 
transition  would  occur  at  lower  values  if  rough¬ 
ness  effects  promoted  transition.  Using  calcu¬ 
lations  at  flight  conditions  corresponding  to  the 
trajectory  point  where  transition  was  assumed 
during  the  TPS  design  it  was  found  that  a  value 
of  Ree/Me  =  250  corresponded  to  a  location  on 
the  windward  centerline  at  80%  of  the  length. 
This  value  was  ultimately  adopted  as  the  X-33 
transition  criterion.  To  insure  that  transition  due 
to  roughness  would  not  occur  earlier  in  the  tra¬ 
jectory,  a  correlation  between  Ree/Me  and 
roughness  height  was  needed  to  determine  a 
maximum  allowable  roughness  for  the  vehicle. 

Berry,  et  al.  [36]  developed  a  successful 
correlation  of  wind-tunnel  data  for  discrete 
roughness  elements  on  the  Shuttle  Orbiter  that 
related  Ree/Me  with  k/8,  where  k  is  the  rough¬ 
ness  height  and  5  is  the  boundary  layer  thick¬ 
ness.  A  similar  approach  was  used  to  obtain  a 
correlation  for  X-33  whereby  raised  square 
patches  of  varying  thickness  were  applied  at 
discrete  surface  locations  and  the  phosphor 
thermography  technique  used  to  determine  the 
“incipient”  and  “effective”  transition  Reynolds 
number.  An  “incipient”  value  corresponded  to 
the  highest  Reynolds  number  that  maintained 
laminar  flow  downstream  of  the  trip  while  an 


“effective”  value  corresponded  to  the  minimum 
Reynolds  number  where  the  transition  front  was 
fixed  at  the  roughness  element.  Tunnel  condi¬ 
tions  were  varied  for  a  given  roughness  height, 
trip  location,  and  angle  of  attack  to  locate  these 
incipient  and  effective  Reynolds  numbers. 
Roughness  elements  were  fabricated  from 
0.0025-inch  Kapton  tape  which  could  be 
stacked  in  multiple  layers  to  provide  thickness 
variation.  The  elements  simulated  the  effect  of 
a  raised  TPS  panel  and  were  oriented  with  the 
diagonal  aligned  to  the  local  flow,  much  like  the 
herringbone  pattern  of  the  metallic  TPS. 
Roughness  elements  were  tested  at  six  locations 
along  the  windward  centerline  and  at  four  loca¬ 
tions  along  the  attachment  line.  Attachment 
lines  on  X-33  are  confined  outboard  near  the 
fuselage  chines  for  the  angles  of  attack  studied 
(20,  30,  and  40-deg)  and  location  of  the  lines 
was  determined  via  oil  flow  images  and  com¬ 
putational  predictions  prior  to  the  subject  tests. 
With  this  arrangement,  a  systematic  testing  of 
roughness  height,  location,  and  angle  of  attack 
was  performed  and  the  incipient  and  effective 
transition  Reynolds  numbers  determined. 
Computing  Ree/Me  and  boundary  layer  thick¬ 
ness  at  each  trip  location  for  all  effective  and 
incipient  pairs  provided  the  necessary  informa¬ 
tion  to  establish  a  correlation  between  roughness 
and  transition.  Figure  17  shows  the  resulting  data 
and  displays  a  generally  well-behaved  trend  in 
terms  of  Ree/Me  vs.  k/8.  The  lines  in  the  figure 
represent  “conservative”  curve  fits  of  the  data 
that  delineate  a  region  of  laminar  flow  below  the 
dashed  line  (the  “incipient”  curve)  and  a  region 
of  turbulent  flow  above  the  solid  line  (the  “ef¬ 
fective”  curve). 
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k/8 


Figure  17.  Experimental  correlation  of  X-33  transi¬ 
tion  with  roughness  from  discrete  trips  (from  Berry,  et 
al.  [37 j). 

The  incipient  curve  shown  in  Figure  17  pro¬ 
vided  a  basis  for  determining  allowable  rough¬ 
ness  heights  on  the  X-33.  Using  this  figure  with 
the  value  of  Ree/Me  =250  (as  determined  for 
X-33)  gives  k/5  =  0.2.  Applying  this  constraint 
(k/5=0.2)  to  the  surface  boundary  layers  calcu¬ 
lated  at  the  transition  trajectory  point  used  in  the 
TPS  design  yielded  the  allowable  heights  over 
the  vehicle  shown  in  Figure  18.  Roughness 
height  less  than  the  allowable  should  be  negligi¬ 
ble  while  larger  values  would  cause  transition 
earlier  than  anticipated.  Comparison  of  the  al¬ 
lowable  heights  to  the  manufacturing  tolerances 
on  X-33  revealed  that  the  roughness  constraints 
were  easily  met  (except  possibly  near  the  fin 
root). 


Figure  18.  Maximum  allowable  roughness  on  full- 
scale  vehicle  based  on  experimental  correlation  of  X- 
33  transition  (from  Thompson,  et  al.  [35]). 

Another  transition  issue  investigated  for 
X-33  was  the  effect  of  a  distributed  roughness 
caused  by  bowing  of  the  metallic  TPS  panels 
expected  in  flight.  This  bowing  is  specific  to 
X-33  due  to  a  rapid  heat  pulse  associated  with 
the  sub-orbital  flight  trajectory;  bowing  would 
not  be  expected  on  an  RLV  during  reentry  from 
LEO.  Figure  19  illustrates  five  configurations  of 
bowed  panel  arrangements  that  were  tested. 
Three  different  bow  heights  were  tested  (0.002, 
0.004,  and  0.008  inch)  with  the  two  smallest 
values  being  geometrically  scaled  from  the  ex¬ 
pected  bowing  for  the  full-scale  vehicle.  The 
wavy  surface  introduced  clear  flow  field  pertur¬ 
bations  observed  in  oil  flow  and  schlieren  im¬ 
ages,  particularly  for  the  largest  two  bow 
heights.  However,  it  was  found  that  all  of  the 
configurations  with  bow  heights  smaller  than 
the  maximum  value  were  less  effective  at  pro¬ 
ducing  transition  than  the  discrete  roughness 
elements.  It  was  concluded  [37]  that  the  bowed 
wall  did  not  cause  any  additional  sensitivity  to 
roughness  beyond  the  correlation  already  devel¬ 
oped  using  discrete  roughness  elements. 
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Figure  19.  Configurations  tested  for  bowed  panel 
roughness  test  (from  Berry,  et  al.  [37]). 

The  simple  Ree/Me  criterion  developed  for 
X-33  enabled  construction  [35]  of  a  semi- 
analytic  model  to  predict  whether  transition  oc¬ 
curs  at  a  station  on  the  windward  centerline  for 
a  given  altitude,  velocity,  and  angle-of-attack 
condition.  This  model  is  based  on  a  database  of 
numerical  simulations  that  cover  the  range  of 
freestream  and  angle  of  attack  values  expected 
during  flight.  For  X-33,  a  domain  in  trajectory 
space  was  selected  where  altitude  ranged  from 
30  to  80  km  for  freestream  velocities  between  1 
and  5  km/s  and  angles  of  attack  between  0  and 
40  deg.  A  total  of  275  flow  solutions  were 
computed  for  discrete  combinations  of  these  pa¬ 
rameters.  Each  calculation  was  performed  using 
the  LATCH  code  to  provide  a  distribution  of 
Ree/Me  along  the  windward  centerline.  Since 
LATCH  requires  edge  properties  from  an  invis- 
cid  flowfield  solution,  it  was  necessary  to  make 
an  approximation  for  the  approach  to  be  practi¬ 
cal  by  using  a  single  inviscid  solution  at  each 
angle  of  attack  for  all  altitudes  and  velocities. 
This  approximation  is  valid  for  all  altitudes  (ex¬ 
cept  when  real  gas  effects  are  important)  since 
the  inviscid  solution  is  independent  of  Reynolds 
number.  It  is  valid  for  velocities  where  the 
Mach  number  is  large  (M>8)  since  the  shock 
shape  and  normalized  pressure  distributions  are 
nearly  invariant.  To  improve  the  accuracy  of 
inviscid  edge  properties  for  the  X-33  trajectory, 


inviscid  equilibrium  air  solutions  were  obtained 
for  Mach  10  cases  at  angles  of  attack  between 
20-  and  40-deg  and  Mach  6  perfect  gas  solutions 
were  used  for  angles  between  0  and  15  deg. 

A  relationship  for  transition  at  a  body  point 
was  derived  by  interpolation  of  the  LATCH  da¬ 
tabase  to  locate  the  altitudes  where  transition 
occurs  for  all  combinations  of  velocity  and  an¬ 
gle  of  attack.  Figure  20  shows  this  relationship 
using  the  X-33  criterion  (Ree/Me=250)  and  a 
point  at  x/L=0.8  on  the  windward  centerline. 
The  result  is  a  nearly  planar  surface  in  three- 
dimensional  trajectory  space  that  separates 
laminar  flow  above  from  turbulent  flow  below. 
In  practice,  this  function  can  be  used  to  identify 
regions  of  transition  by  plotting  a  trajectory  in 
this  space  and  locating  the  point  where  the  path 
intersects  the  surface  as  shown  in  the  figure.  In 
addition,  the  sensitivity  of  transition  onset  to  the 
value  of  Ree/Me  and  the  movement  of  transition 
along  the  vehicle  can  be  studied  by  varying  the 
parameters  used  for  interpolation  in  the  data¬ 
base.  A  computer  routine  to  evaluate  the  transi¬ 
tion  model  for  X-33  was  developed  that  pro¬ 
vided  invaluable  information  to  the  vehicle  de¬ 
signers  leading  to  tailoring  of  trajectories  in  or¬ 
der  to  remain  within  TPS  design  limits. 


Figure  20.  Transition  criterion  in  trajectory  space 
(from  Thompson,  et  al.  [35]). 
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5  Summary 

Development  of  a  sub-scale  flight  vehicle 
(X-33)  to  demonstrate  technologies  for  a  new 
reusable  launch  system  in  the  United  States  has 
ensued  over  the  past  five  years.  The  new  tech¬ 
nologies  and  the  rapid  program  pace  have  pre¬ 
sented  many  challenges  to  industry  and  gov¬ 
ernment  teams.  In  this  paper,  an  overview  of 
several  major  studies  conducted  by  the  Aero- 
thermodynamics  Branch  at  NASA  LaRC  to 
support  development  of  the  LMSW  X-33  vehi¬ 
cle  is  presented.  These  studies  provided  aero¬ 
dynamic  and  aeroheating  screening,  optimiza¬ 
tion,  and  benchmarking  of  the  X-33  and  utilized 
the  unique  capabilities  offered  by  the  LaRC 
Aerothermodynamics  Facilities  Complex  cou¬ 
pled  with  the  application  of  computational  tools 
(both  detailed  CFD  and  engineering  level). 

Extensive  wind  tunnel  testing  in  the  LaRC 
Mach  6  and  10  air  tunnels  and  Mach  6  CF4  tun¬ 
nel  contributed  significantly  to  the  establish¬ 
ment  of  the  flight  databook  for  the  hypersonic 
longitudinal  and  lateral  aerodynamic  character¬ 
istics,  control  surface  effectiveness,  and  reaction 
control  system  effects.  Measurements  from 
these  tests  showed  the  baseline  vehicle  (zero 
control  deflections)  to  be  marginally  stable  in 
longitudinal  pitch  and  laterally  unstable  at  the 
high  trim  angles  of  attack  (45  to  50  deg). 
Measurement  and  prediction  revealed  the  vehi¬ 
cle  became  less  stable  as  Mach  number  de¬ 
creased  from  10  to  6  and  the  source  of  this  trend 
was  isolated  to  the  vehicle  fuselage.  While  this 
Mach  number  effect  was  qualitatively  predicted 
by  computation  and  the  measured  and  predicted 
pitching  moment  were  in  good  agreement  at 
Mach  10,  the  prediction  at  Mach  6  fell  outside 
the  experimental  uncertainty. 

Wind  tunnel  tests  with  deflected  body  flap 
surfaces  showed  the  vehicle  to  have  more  than 
sufficient  authority  for  pitch  control,  which 
mitigated  concerns  about  the  decreasing  stabil¬ 
ity  with  Mach  number.  Body  flap  effectiveness 
suffered  degradation  with  increasing  deflection 
angle  and  increasing  angle  of  attack  due  to  a 
shock  interaction  and  impingent  of  a  reflected 
expansion  fan  on  the  flap  surface.  Predictions 


with  CFD  codes  significantly  enhanced  the  un¬ 
derstanding  of  the  complex  flowfield  around  the 
deflected  flap  and  were  in  good  comparison 
with  the  measured  pitching  moments.  Extensive 
wind  tunnel  testing  on  an  X-33  model  with  re¬ 
action  control  jets  was  performed  to  establish  a 
database  for  the  aerodynamic  effect  of  jet  inter¬ 
action  with  the  flowfield  and  vehicle  surface. 
Side  firing  jets  were  moved  to  the  base  region  as 
a  result  of  these  tests. 

Global  heat  transfer  distributions  measured 
using  a  thermographic  phosphor  technique  in 
the  Mach  6  air  tunnel  were  in  good  agreement 
with  predictions  from  the  CFD  codes  used  to 
define  the  aerothermal  environment  in  flight. 
Detailed  studies  of  the  aerodynamic  heating  to 
the  deflected  flap  surfaces,  aerospike  engine 
nozzle,  and  fin  leading  edge  were  performed 
with  the  phosphor  technique  to  provide  infor¬ 
mation  about  these  critical  regions  to  the  vehicle 
designers.  Predicted  engine  nozzle  heating  due 
to  an  impingement  from  the  wind-side  flow 
field  at  flight  conditions  compared  reasonably 
well  with  wind  tunnel  measurements.  These 
values  of  predicted  heating  and  temperatures  in 
flight  contributed  significantly  to  the  decision  to 
employ  active  cooling  to  the  engine  nozzles  us¬ 
ing  hydrogen  fuel  during  unpowered  descent. 

Global  aeroheating  tests  at  Mach  6  over  a 
range  of  Reynolds  number  and  incidence  for  a 
“smooth  body”  with  and  without  discrete 
roughness  trips  along  the  windward  centerline 
and  outboard  attachment  lines  were  correlated  to 
give  a  relation  between  the  boundary  layer  tran¬ 
sition  parameter  (Ree/Me)  and  the  roughness 
height.  Applying  this  roughness  correlation  at 
flight  conditions  produced  a  surface  map  of  the 
maximum  allowable  roughness  heights  for  the 
full-scale  vehicle.  Fabrication  and  testing  of 
thermographic  phosphor  models  with  “wavy- 
wall”  surfaces  to  simulate  bowing  of  metallic 
themal  protection  tiles  in  flight  was  performed. 
These  tests  revealed  that  the  discrete  trips  were 
more  efficient  generators  of  boundary  layer 
transition  than  the  wavy  wall.  The  transition 
parameter  established  in  the  experimental  work 
was  used  to  derive  a  semi-analytic  model  capa¬ 
ble  of  predicting  transition  for  X-33  at  points  on 
the  windward  fuselage  over  a  wide  range  of  al- 
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titude,  velocity,  and  angle  of  attack.  Vehicle  de¬ 
signers  used  this  model  to  tailor  trajectories  so 
not  to  exceed  the  TPS  design  limits. 

This  review  of  studies  performed  at  NASA 
LaRC  to  support  the  full  breadth  of  hypersonic 
aerodynamic  and  aerothermodynamic  issues  for 
X-33  demonstrates  the  advantages  and  credibil¬ 
ity  afforded  by  synergism  between  experimental 
and  computational  capabilities.  The  ability  to 
analyze  critical  aerodynamic  and  aerothermody¬ 
namic  design  issues  is  greatly  enhanced  by  the 
synthesis  of  these  disciplines. 
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